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Abstract: Feynman diagrams with two real partons contributing to the 
next-to-leading-order singlet gluon-quark DGLAP kernel are analysed. The 

infra-red singularities of unintegrated distributions are examined numeri- 

^ r~| cally. The analytical formulae are also given in some cases. The role of 

the colour coherence effects is found to be crucial for cancellations of the 
I double- and single-logarithmic infra-red singularities. 

Qh 



> 

X 



Submitted To Acta Physica Polonica B 
PACS numbers: 12.38.-t, 12.38.Bx, 12.38.Cy 



> 

o 

IFJPAN-IV-2011-5 

o 



* This work is supported by the Polish Ministry of Science and Higher Education grant 

No. 1289/B/H03/2009/37. 

Presented by Magdalena Slawinska at the Cracow Epiphany Conference 2011 - on 
the First Year of the LHC, January 10-12, 2011 



(1) 



2 



IFJPAN-IV-2011-5 PRINTED ON October 5, 2011 



1. Motivation 

The study presented here is a part of the development of the fuUy ex- 
clusive next-to-leading order (NLO) Parton Shower Monte Carlo (MC) for 
precision QCD predictions for the LHC experiments, see (l]-[3]. DGLAP [4] 
evolution of parton distribution functions (PDFs) is modelled in the Monte 
Carlo within the unintegrated phase space. A methodology based on the 
collinear factorisation theorems in physical gauge based on refs. [5] and [g] 
IS usedQ The MC program will simulate exactly NLO DGLAP evolution 
of PDFs by itself, as opposed to using pretabulated PDFs, provided by 
the non-MC programs like QCDNUM |7|. For the construction of such a 
new NLO parton shower MC program a new exclusive (fully unintegrated) 
NLO evolution kernels are required in order to impose NLO corrections, by 
means of reweighting the LO distribution, as outlined in ref. [3]. In this 
method LO MC parton shower has to be reconstructed from the scratch, 
contrary to methodology of ref. [IT], where at the price of MC weights being 
non-positive, one is able to use standard LO parton shower MC. 

Very schematically, the corresponding MC weight reads: 

, , exact NLO diagram distribution 

MC weight = , — : . (1) 

crude LO distribution 

The potential problem is that a single Feynman diagram, or small subset of 
diagrams, entering the exclusive NLO kernel (and the MC weight) generally 
is not gauge-invariant and may feature uncancelled soft singularities. The 
Monte Carlo weights may then explode, unless the crude distribution of the 
LO MC already reproduces exactly soft singularities of NLO diagrams. It 
is therefore very important to understand in fine detail the structure of soft 
and collinear singularities in exclusive kernels, as implemented in the LO 
MC and also including complete NLO corrections. 

For the purpose of the MC we are going to analyse the example dia- 
grams and their groups one by one, gaining the detailed knowledge about 
the structure of collinear and soft singularities of each Feynman diagram 
contributing to the NLO kernel and the interplay between diagrams. We 
shall exploit tools and methods of the graphical analysis of the infra-red sin- 
gularities which were already used for the non-singlet diagrams in ref. }12\. 
Here, we will extend this study to a gauge invariant subset of two-real singlet 
diagrams contributing to the Pgq NLO DGLAP kernel. Let us stress that 
the cancellations discussed in the following are not of the usual KLN fls] na- 
ture, i.e. between the real and virtual Feynman diagrams, but rather among 
the real diagrams alone, and are governed by the spin and colour quantum 



^ In the complementary approach of refs. [8| |l0| soft singularities are resummed first 
and collinear resummation is added next. 
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numbers. The contributions of the diagrams to the standard DGLAP (in- 
clusive) kernels analysed in the following have been already defined and used 



in refs. |14;| and 15 . Generally, we shall examine the structure of the soft 
and collinear singularities of the unintegrated distributions related to these 
NLO DGLAP evolution kernels. 



2. Singlet diagrams considered 

The singlet diagrams considered in this contribution originate from the 
LO amplitude for splitting a gluon into a quark (antiquark) 



by means of adding the NLO corrections from the emission of an additional 
gluon: 




Fig. 1. Singlet gluon-quark diagrams; "1" denotes a quark and "2" - a gluon. 



Feynman diagrams contributing to the NLO kernel result from squaring 
the above sum of amplitudes and are displayed in Fig. [l} It is worth noting, 
that the ladder diagrams (the first and the third in the upper row in Fig. [T]) 
enter the NLO kernel supplemented with the so called collinear countert- 
erms that subtract off the leading-order contributions. In this contribution, 
however, the counterterms will be included only at the end of the analysis, 
and before that the leading-order singularities will be visible. 
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Fig. 2. Amplitude-squared diagrams ^ Cp 



We will adopt the same approach as in ref. [12] and analyse all distri- 
butions in the logarithmic Sudakov variables (ln(ai/a2), ln(ai/a2)), where 
ai come from the Sudakov parametrisation of four-momenta of the emitted 
particles: ki = aip + a~n + ki± and ai are angular (rapidity-related) vari- 
ables, Oj = i = 1,2. All contributions are normalised to the eikonal 
phase space: 

dai da2 dai da2 
d^- = , (2) 

ai a2 ai a2 

with the angles integrated over and terms of order 0(e) neglected. More- 
over, we will ensure that at least one emission is hard by constraining 
ai + 02 = 1 — X > 0. Similarly, the maximal angle is fixed to an arbitrary 
parameter. Hence, the ratios ai/a2 and ai/a2 will measure the relative 
hardness and angles of the two partons. 

We will explore the soft limit of the diagrams in Fig. [T| namely the limit 
where both |kj_|_| — )• and — )• (for a given i), but Oj remains finite. In 
the logarithmic Sudakov variables the singularities will appear on the plots 
as one- or two-dimensional infinite structures. 



3. Results 

Let us first consider the Cp class of diagrams, corresponding to emission 
of a gluon from a quark. The bremsstrahlung type diagram is displayed in 
Fig. ^(left). It has a doubly- logarithmic singularity visible as the infinite 

trapezoidal plateau, bordered by the lines ai = 02 and ^ = ^ (the line of 

equal "lightcone minus variables" a~). The second diagram, representing 
the amplitude-squared of the emission of a gluon from the emitted quark 
(Fig. [2| middle), features a collinear singularity manifesting itself as the 
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infinite ridge along the line of equal angles. This singularity, however, is 
not related to the soft limit (it is compensated by the virtual diagram) and 
will not be considered here. This diagram has also a doubly-logarithmic 
singularity in the form of a triangular plateau, bordered by the lines of 
equal angles and "minus variables" . The sum of the two (the rightmost plot 
in Fig. [2]) features two equal-height plateaux with the canyon at the line of 
equal minus variables. 

The most singular terms from the distributions of the diagrams are nec- 
essarily proportional to the products of the leading-order DGLAP kernels: 



Q'^(ai,a2) 



(3) 



and 



P,giZl)P,,iZ2) = CP 



2 af + (1 — aiY -I- (1 — «! 



'1 — ai)a2 



(4) 



where zi = 1 — ai and Z2 



1— ai 



Similarly: 



2"l(^l+(l_ 
p 9 — 



a? g'*(ai,a2)' 



(5) 



P,g{zl}P,,{Z2) - Cp - 

where zi = 1 — x and Z2 = jz^- 

In eqs. ^ and ([s]) we also used (?^(ai, 02) = "^1+ a\+2aia2 cos (/> 
(proportional to the denominator of the propagator of the most virtual 
quark) and a? = a\ + a^ — 2aia2 coscp (proportional to the invariant mass 
of the emitted quark and gluon). 

The "canyon" structure in the plot, being the remaining singly logarith- 
mic singularity, however, spoils the soft limit regardless of the counterterm 
employed. If we add now the interference diagramj^ as shown in Fig. [sj the 
canyon gets removed. 

What remains is the uniform plateau bordered by ai = 02 coUinear 
singularity. The "minus variable" ordering, preferred by each diagram sep- 
arately, turns out to be irrelevant for the sum of diagrams! 



This diagram's colour coefficient is equal to Cp — CaCf/2- In this analysis we add 
only its part ~ Cp, adding the other one ~ —CaCp/2 to the CaCf diagrams. 
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Fig. 3. The infra-red cancellations among the diagrams ^ Cp. 



Prom the above formulae it follows that the quadratic plateau represents 
the leading-order contribution. In the NLO kernel it is removed by the 
counterterm of the factorisation procedure proportional to 



2{x' + {1 - a,r) (aj + {I 



ai 



1 — «! 



1 — «! 



(7) 

The doubly logarithmic singularity of the counterterm depends on the de- 
tails of factorisation procedure in usej^ 

The CaCf subset consists of diagrams that correspond to the emission 
of a gluon from the incoming gluon. They include only one amplitude- 
squared diagram. Displayed in Fig. |4j left, it has a doubly-logarithmic 
singularity - the plateau stretching in the two regions: where the angle of 
the emitted gluon is larger than the quark's and vice versa. The diagram 
contributes in the region of phase space where both emissions are ordered in 
the minus variable. After adding the interference diagram, the boundaries 
of the resulting plateau are corrected and the sum contributes in the region 
of the phase-space, where the angle of the quark is larger than the angle of 
the gluon. 

The collinear counterterm is given by: 

(1 - 02)^ (1 - 02) 1-02 



(8) 



The counterterm in eq. ([7| has the additional theta function related to ordering of 
the emissions (not shown explicitly), defining the boundaries of the LO plateaux. In 
the following we will use the ordering in the angular variables a^. 
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Fig. 4. The cancellation of doubly-logarithmic singularities among diagrams ^ 
CfCa- 




Fig. 5. The cancellation of singly-logarithmic singularities among interference dia- 
grams ~ CfCa- 

While the Cp counterterm of eq. Q features the oi < 02 ordering, the 
ordering in the CaCf counterterm of eq. ([s]) is the opposite (02 < ai) due 
to a gluon being emitted before a quark. 

The remaining CaCf interferences feature single-log singularities, seen 
as infinite canyons/ridges along the hne of equal minus- variables in Fig. [5] 
that cancell out when added. What remains is the little hill in the central 
region, which leads to a finite contribution. 

The sum of all singlet diagrams discussed in this contribution is pre- 
sented in Fig. |6j The left-hand-side plot in this figure shows two leading- 
order plateaux separated by the line of equal angles ai = 02- The line 
represents a collinear singularity and comes from the diagram, in which the 
additional gluon is emitted from the emitted quark. The plateau on the left 
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Fig. 6. All singlet gluon-quark contributions to the NLO kernel added together 
(left) and with counterterms subtracted (right). 

(in brown) corresponds to the topologies in which a soft gluon is emitted 
from a quark. The right-hand-side plateau (navy-blue) represents contri- 
butions with a soft gluon emitted from the incoming gluon. The relative 
height of both plateaux is equal to Cp/CpCA, as expected. 

In the right plot the same sum is presented, but with the leading or- 
der singularities cancelled out by the factorisation counterterms of eqs. ([7]) 
and ^ on the left- and right- hand-side of this plot, respectively. The plot 
features the collinear singularity only. 

4. Conclusions 

We conclude that the restoration of gauge invariance (colour coherence) 
is crucial in cancelling infra-red singularities. We understand the soft limits 
of NLO exclusive kernels, observe and explain the cancellations of double- 
and single-logarithmic soft singularities. The angular ordering is the pre- 
ferred parametrisation of the phase space in view of the soft singularity 
structure of the distributions from gauge-invariant subset of diagrams con- 
tributing to NLO evolution kernels in the exclusive (unintegrated) form. 

REFERENCES 

[1] S. Jadach, M. Skrzypek, A. Kusina, and M. Slawinska, Pes' RADCOR2009 
(2010) 069, 1002.0010 

[2] S. Jadach, A. Kusina, M. Skrzypek, and M. Slawinska, | 1102 . 508 31 

[3] S. Jadach, A. Kusina, W. Placzek, M. Skrzypek, and M. Slawinska, 
11103.50151 



IFJPAN-IV-2011-5 printed on October 5, 2011 



9 



4] L.N. Lipatov, Sov. J. Nucl. Phys. 20 (1975) 95; 

V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; 
G. Altarelli and G. Parisi, Nucl. Phys. 126 (1977) 298; 
Yu. L. Dokshitzer, Sov. Phys. JETP AQ (1977) 64. 

5] R. K. Ellis, H. Georgi, M. Machacek, H. D. Politzer, and G. G. Ross, Phys. 
Lett. B78 (1978) 281. 

G. Curd, W. Furmanski, and R. Petronzio, Nucl. Phys. B175 (1980) 27. 
M. Botje, Comput. Phys. Commun. 182 (2011) 490-532, _1005 . 1481, 



B. F. L. Ward, Annals Phys. 323 (2008) 2147-2171, 0707 . 3424 



S. Joseph, S. Majhi, B. F. L. Ward, and S. A. Yost, Phys. Lett. B685 (2010) 
283-292, 0906.0788, 

S. Jose ph, S. Ma jhi, B. F. L. Ward, and S. A. Yost, Phys. Rev. D81 (2010) 
076008, fl00l7l 434, 

S. Frixione and B. R. Webber, JHEP 06 (2002) 029, |hep-ph /0204244 

M. Slawinska and A. Kusina, Acta Phys. Polon. B40 (2009) 2097-2108, 
0905.1403 

T. Kinoshita, J. Math. Phys. 3 (1962) 650; 

T.D. Lee and M. Nauenberg, Phys. Rev. 133B (1964) 1549. 

W. Furmanski and R. Petronzio, Phys. Lett. B97 (1980) 437. 

R. K. Ellis and W. Vogelsang, hep-ph/9602356 



